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ABSTRACT 



We analyze MOA-2010-BLG-311, a high magnification (A max > 600) microlensing 
event with complete data coverage over the peak, making it very sensitive to planetary 
signals. We fit this event with both a point lens and a 2-body lens model and find that 
the 2-body lens model is a better fit but with only Ax 2 ~ 140. The preferred mass 
ratio between the lens star and its companion is q = lO -3 ' 71 * 10 ' 1 , placing the candidate 
companion in the planetary regime. Despite the formal significance of the planet, we 
show that because of systematics in the data the evidence for a planetary companion 
to the lens is too tenuous to claim a secure detection. When combined with analyses 
of other high-magnification events, this event helps empirically define the threshold for 
reliable planet detection in high- magnification events, which remains an open question. 

Subject headings: Galaxy: bulge — gravitational lensing: micro — planetary systems: 
detection 



1. Introduction 

High-magnification events, events in which the maximum magnification of the source, A mSL ^, is 
greater than 100, have been a major focus of microlensing observations and analysis. Because the 
impact parameter between the source and the lens is very small in such cases, uq ~ l/Anax, it is 
likely to probe a central caustic produced by a planetary companion to the lens star. Furthermore, 
such events can often be predicted in advance of the peak, allowing intensive observations of the 
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event at the time when it is most sensitive to planets. Consequently, a substantial amount of effort 
has been put into identifying, observing, and analyzing such events. 

Observed high-magnification events are classified into two groups for further analysis: events 
with signals obvious to the eye and events without. Only events in the first category are system- 
atically fit with 2(or more)-body models. The other events are only analyzed to determine their 
detection efficiencies. As a result, no planets have been found at or cl ose to the de t ection thresh- 
old, and furthermore this detection threshold is not well understood^. iGould et al. (201C) suggest 



a detection threshold in the range of A% 2 = 350-700 is required to both detect the signal and 
constrain it to be planetary, but they note that the exact value is unknown. With the advent of 
second-generation microlensing surveys, which will be able to detect planets as part of a controlled 
experiment with a fixed observing cadence, it is important to study the reliability of signals close to 
the detection threshold, since a systematic analysis of all events in such a survey will yield signals 
of all magnitudes, some of which will be real and some of which will be spurious. 

In this paper, we present the analysis of a high-magnification microlensing event, MOA-2010- 
BLG-311, which has a planetary signal slightly too small to claim as a detection. We summarize 
the data properties in Section [2] and present the color-magnitude diagram (CMD) in Section [3l In 
Section [U we fit the light curve with both point lens and 2-body models. We then discuss why 
a planetary detection cannot be claimed in Section 14.41 We calculate the Einstein ring size and 
relative proper motion in Section [5] and discuss the possibility that the lens is a member of the 
cluster NGC 6553 in Section [6l We give our conclusions in Section [71 

2. Data 
2.1. Observations 

On 2010 June 15 (FDD' 53 62.967 = HJD-2450000), the Micr olensing Observations in As- 



trophysics (MOA) collaboration ([Bond et al.ll200ll ; ISumi et al.ll201ll ) detected a new microlensing 



event MOA-2010-BLG-310 at (R.A., decl.) = (18 h 08 m 49. s 98, -25°57'04."27) (J2000.0), (1, b) = 
(5.17, -2.96), along our line of sight toward the Galactic Bulge. MOA announced the event through 
its email alert system and made the data available in real-time. Within a day, this event was 
identified as likely to reach high magnification. Because of MOA's real-time alert system, the event 
was identified sufficiently far in advance to allow intensive follow up observations over the peak. 

The observational data were acquired from multiple observatories, including members of the 
MOA, OGLE, /iFUN, PLANET, RoboNet, and MiNDSTEp collaborations. In total, sixteen ob- 
servatories monitored the event for more than one night, and thus their data were used in the 
following analysis. Among these, there is the MOA survey telescope (1.8m, MOA-RecH, New 



The need for a well-defined detection threshold is also discussed in lYee et al-llioH 
2 This custom filter has a similar spectral response to i?-band 



- 6 - 



Zealand) and the B&C telescope (60cm, V, I, New Zealand); eight of the observatories are from 
/iFUN: Auckland (AO, 0.4m, Wratten #12, New Zealand), Bronberg (0.36m, unfiltered, South 
Africa), CTIO SMARTS (1.3m, V, I, H, Chile), Farm Cove (FCO, 0.36m, unfiltered, New Zealand), 
Kumeu (0.36m, unfiltered, New Zealand), Molehill (MAO, 0.3m, unfiltered, New Zealand), Vintage 
Lane (VLO, 0.4m, unfiltered, New Zealand), and Wise (0.46m, unfiltered, Israel); three are from 
PLANET: Kuiper telescope on Mt. Bigelow (1.55m, I, Arizona), Canopus (1.0m, /, Australia), 
and SAAO (1.0m, I, South Africa); one is from RoboNet: Liverpool (2.0m, I, Canaries); and one 
is from MiNDSTEp: La Silla (1.5m, I, Chile). The event also fell in the footprint of the OGLE IV 
survey (1.3m, /, Chile), which was in the commissioning phase in 2010. The observatory and filter 
information is summarized in Table [TJ 

In particular, observations from the MOA survey telescope, MOA B&C, PLANET Canopus, 
/iFUN Bronberg, and /iFUN VLO provided nearly complete coverage over the event peak between 
HJD' = 5365.0 and HJD' = 5365.4. 



2.2. Data Reduction 



The MOA and B&C data were reduced with the standard MOA pipeline (jBond et al 



2001 



The d ata from the /iFUN observatories were reduced using the standard DoPhot reduction (jSchechter et al 
19931 ). with th e exception of Bronberg and VLO data, which were reduced using difference image 

dl boool 



analysis (DIA; lAlar 



Wozniak 2000). Data from the PLANET and RoboNet collaborations 



were reduced using pySIS2 pipeline (jBramichl 120081 ; lAlbrow et al J 120091 ). Data from MiNDSTEp 
were also initially reduce d using the py SIS2 pipeline. The OGLE data were reduced using the 
standard OGLE pipeline (judalskil I2OO3I ) . Both the MOA and MiNDSTEp/La Silla data were re- 
duced in real-time, and as such the initial reductions were sub-optimal. In fact, the original MOA 
data over the peak were unusable because they were corrupted. After the initial analysis, both the 
MiNDSTEp and MOA data were rereduced using optimized parameters. 



2.3. Error Renormalization 

With microlensing data, it is frequently the case that the errors calculated for each data point 
by the photometry packages underestimate the true errors. Additionally, there can be outliers in 
the data that are clearly seen to be spurious by comparison to other data taken simultaneously from 
a different site. Simply taking the error bars at face value would lead to biases in the modeling. 
Because the level of systematics varies between different data sets, underestimated error bars can 
give undue weight to a particular set of data. Additionally, if the errors are underestimated, the 
relative A% 2 between two models will be overestimated, making the constraints seem stronger than 
they actually are. 

To first order, we can compensate for the underestimated error bars by rescaling them by a 
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single factor, k. The rescaling factor is chosen for each data set, i, such that ki = (x^/iVj) 1 / 2 , pro- 
ducing a x 2 P er degree of freedom of Xd Q { = 1- We will refer to this simple scheme for renormalizing 
the errors as "1-parameter errors". 

Alternatively, we can use a more complex method t o renormalize the e rrors, which we will 



call "2 -parameter errors". This method was also used in iMiyake et al.l ( 2012l ) and iBachelet et al 



(|2012l ). To renormalize the errors, we rank order the data by magnification and calculate two 



factors, k and e m i n , such that 



cTj — ki^al rigj + e^ ini , (1) 

where o" or ig is the original error bar and a is the new error bar and the calculations are done in 
magnitudes. The index i refers to a particular data set, and j refers to a particular point within that 
data set. The additional term, e m i n , enforces a minimum uncertainty in magnitudes, because at 
high magnification, the flux is large, so the formal errors on the measurement can be unrealistically 
small. The error factors, k and e m i n , are chosen s o that y\ r for points sorted by magnification 



increases in a uniform, linear fashion and Xdof = 1 ( Yee et al.l 120121 ). Not all data sets will require 
an e min term; it is only necessary in cases for which Xdd ^ as a break because the formal errors 
are too small when the event is bright. Note that the e m ; n factor will be primarily affected by the 
points taken when the event is bright, whereas k is affected by all points, so if there are many more 
points at the baseline of the event, these will dominate the calculation of k. 

We also remove outliers in the data. We begin by eliminating any points taken at airmass > 3 
or during twilight. Additional outliers are defined as any point more than Xa from the expected 
value, where X is determined by the number of data points, such that fewer than one point is 
expected to be more than Xa from the expected value assuming a Gaussian error distribution. 
The normal procedure is to compare the data to the "expected value" from a point lens fit. We do 
this for data in the wings of the event, i(HJD') < 5363 or t(HJD') > 5367, when we do not expect 
to see any real signals. However, the peak of the event, 5365 < i(HJD') < 5367, is when we would 
expect to see a signal from a planet if one exists. A planetary signal would necessarily deviate from 
the expectation for a point lens, so in this region instead of comparing to a point lens model, we 
use the following procedure to identify outliers: 



1. For each point, we determine whether there are points from any other data set within 0.01 
days. If there are more than two points from a given comparison dataset in this range, we 
keep only the point(s) immediately before and/or after the time of the point in question (i.e., 
a maximum of 2 points). The point is not compared to other points from the same data set. 

2. If there are matches to at least two other data sets, we proceed to ([3]) to determine whether 
or not the point is an outlier. Otherwise, we treat the point as good. 

3. We then determine the mean of the collected points, fx, including the point in question, by 
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maximizing a likelihood function for data with outliers (jSivia &: Skilling 1120101 ): 

N 

(2) 

where Rj = (// — Xj)/aj and Xj is the datum and Uj is its renormalized error bar. If the flux 
is changing too rapidly for the points to be described well by a mean, we fit a line to the 
data. 

4. We compare the point to the predicted value, using the likelihood function to see if it is an 
outlier: 

(a) If there is only one likelihood maximum, we calculate Rj. If Rj > X, the point is rejected 
as an outlier. 

(b) If there is more than one maximum, and the point in question falls in the range spanned 
by the maxima, we assume the point is good. If it falls outside the range, we calculate -Re- 
using the nearest maximum to determine ji. If Rj > X, we flag the point as an outlier. 



This procedure is more complicated than usual, but because we compare the points only to other 
data, rather than some unknown model, it provides an objective means to determine whether a 
point is an outlier without destroying real signals corroborated by other data. We also visually 
inspect each set of points to confirm that the algorithm works as expected. 

We iteratively determine the error renormalization factors and outliers. We fit the data to a 
point lens and find the error renormalization factors. We remove the outliers according to the above 
procedure based on the renormalized errors, refit, and recalculate the error renormalization factors. 
We repeat this process until no further outliers are found. For both 1-parameter and 2-parameter 
errors, Table [1] lists the error normalization factors and number of points for each observatory that 
survived the rejection process. 



3. Color-Magnitude Diagram 

To determine the intrinsic source color, we construct a color-magnitude diagram (CMD) of the 
field of view containing the lensing event (Figure [1]) based on V- and /-band images from CTIO 
SMARTS ANDICAM camera. The field stars in the CMD are determined from three F-band 
images and multiple /-band images. Four l/-band images were taken; however, only three of the 
images are of sufficient quality to contribute to the CMD. We visually checked each of the three 
images to make sure that there were no obvious defects such as cosmic-ray events in the images. 

From the fit to the light curve, we find the instrumental magnitude of the source is / = 19.0 
with an uncertainty of 0.05 mag due to differences between the planet and point lens models as 
well as the error parameterizations. In the top right-hand corner of Figure [IJ the red dot at 
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(V — I, I) c \ = (0.40, 15.5) marks the centroi d of the red clump . The intrinsic col or and magnitude 
of the clump are (V-/, J) d)0 = (1-06, 14.3) (jBensbv et alJboill : E^ataf et alJboij ). Using the offset 
between the intrinsic magnitude of the clump and the observed, instrumental magnitude, we can 
then calibrate the magnitude of the source to find Is,o = 17.8 ± 0.1. 

The col or of the source is normally estimated from V- and /-band images using the standard 
technique in lYoo et al.l (|2004l ). However, with only one highly-magnified U-band image, this method 
is unreliable, so we use an alternative technique to determine the instrumental (V — I) color by 
converting from the instrumental (/ — H) color. Using the simultaneous CTIO I- and i7-band 
observations, we measure the instrumental (I — H) color of the source by linear regression of H on 
I flux at various magnifications during the event. We then construct a VIH instrumental color- 
color diagram from stars in the field (bottom panel of Figured]). The stars are chosen to be all 
stars seen in all three bands with instrumental magnitude brighter than Hqtio = 19-0 (note that 
the field of view for fZ-band observations is 2.4' x 2.4' compared to 6' x 6' for the optical bands). 
The field stars form a well defined track, which enables us to estimate the (V — I) source color 
from the observed (I — H) source color. This yields (V — I)o = 0.75 ± 0.05. Note that this method 
would not work for red star s, (V — I)q > 1.3, bec ause for these red stars, the VIH relation differs 
between giants and dwarfs (jBessell fe Brettlll988l ). However, the observed color is we ll blueward of 
this bifurcation. There is also a spectrum of the source taken at HJD' = 5365.001 (jBensbv et al. 
). The "spectroscopic" (V — J)o reported in that work is 0.77, in good agreement with the 



2011 



value calculated here. 



4. Modeling 



4.1. Preamble 



We know that there are correlated errors in microlensing data, but there has not been a 
systematic investigation of this in the microlensing literature. However, correlated errors have a 
major impact on our ability to determine whether or not the planetary signal in this event is real. 

Correlated errors (red noise) are generally thought of as reducing the sensitivity to signals, 
because successive points are not independent, giving related information. But in fact, sharp, 
short timescale signals are not degraded by correlated noise and may still be robustly detected. 
Consider the case of large-amplitude red noise superimposed on small-amplitude white noise (i.e., 
statistically uncorrelated noise). If we set the error bars by the large-amplitude deviations, which 
is correct for long timescales for which the data are uncorrelated, a short timescale jump may be 
considered statistically insignificant. However, if the timescale of the signal is much shorter than 
the correlation length of the red noise, the significance of the signal should be judged against the 
white noise, since on that timescale, the red noise will only contribute a constant offset. 

In principle, the noise properties of the data are fully described by a covariance matrix of all 
data points, but we are unable to calculate such a matrix. Instead, we have two different error 
parameterizations. The k factor is calculated primarily based on baseline data for which statistical 
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errors dominate. In contrast, the e m ; n factor is heavily influenced by points at the peak of the 
event when correlated noise is important. Thus, 2-parameter errors better reflect the red noise, 
whereas 1-parameter errors better reflect the white noise. Because the source in this event crosses 
the position of the lens, finite source effects are pronounced and we expect that any planetary 
signals due to small caustics will occur at the limb-crossing times (tii m b = io i *E \/ P 2 — u q) during 
the peak of the event and last for only ~ 15 minutes (the duration of the signal is equal to £e 
times the caustic size w). In this case, the signal lasts for 10 minutes, which is likely much shorter 
than the timescale for the correlated noise. For typical events, we sample the light curve with a 
frequency of ~ 15 minutes and observe correlations on timescales longer than that, e.g., hours. 
Since the timescale of the signal is less than the timescale of the correlated noise, we begin by 
analyzing the light curve using 1-parameter errors. We will then discuss how the situation changes 
for 2-parameter errors. 

4.2. Point Lens 



A casual inspection of the light curve does not show any deviations from a point lens, so we 
begin by fitting a point lens model to the data. A point lens model is characterized by three basic 
parameters: the time of the peak, t$, the impact parameter between the source and the lens stars, 
uo, and the Einstein timescale, Since uq is small, the finite size of the source can be important. 
To include this effect, we introduce the source size in Einstein radii, p, as a parameter in the 
model. Additionally, we include limb-darkening of the source. The tempe rature of this s lightl y 
evolve d source was determined from the spectrum to be T e fj ~ 5460 K by iBensbv et al.l (|201ll ) . 
Using IClaretl (|2000l ) we found the limb-darkening coefficients to be uy = 0.7086, uj = 0.5470, and 
uh = 0.3624, assuming a microturbulent velocity = 1 k m s -1 , log q = 4.0, solar metallicity, and 



T e ff = 5500K, which is the closest grid point given th e IBensbv et a 



then convert uy, ui, and uh to the form introduced by lAlbrow et al. 



(2011) measurements. We 



2u 



(3) 



to obtain Ty = 0.62, Tj = 0.45, and Tjj = 0.28. Because the various data sets are not on a common 
flux scale, there are also two flux parameters for each data set, /s,i and /b,i, such that 



/mod,i = /s,t^(*) + Ib, 



(4) 



where A(t) is the predicted magnification of the model at time, t. The source flux is given by /g^, 
and /b,« is the flux of all other stars, including any light from the lens, blended into the PSF (i.e., 
the "blend"). 

The best-fit point lens model and the uncertainties in the parameters are given in Table [2j 
This model is shown in Figure El and the residuals to this fit are shown in Figure [3l These exhibit 
no obvious deviations. These fits confirm that finite source effects are important, since p is well 
measured and larger than the impact parameter, uq. 
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We also fit a point lens model that includes the microlens parallax effect, which arises either 
from the orbital motion of the Earth during the event or from the difference in sightlines from two 
or more observatories separated on the surface of the Earth. Microlens parallax enters as a vector 
quantity: 7Te = (tte,Nj tte,e)- The addition of parallax can break the degeneracy between solutions 
with uq > and u < 0, so we fit both cases. Parallax does improve the fit beyond what is expected 
simply from adding 2 more free parameters and shows a preference for no > 0. However, we shall 
see in the next section that a planetary model without parallax produces an even better fit and 
adding parallax in addition to the planet gives only a small additional improvement. 

4.3. 2-body Models 

We search for 2-body models over a broad range of mass ratios, from q = 10~ 6 to q = 10 _1 . 
For each value of q, we chose a range for the projected separation between the two bodies, s, for 
which the resulting caustic is smaller than p and s < 1.0. For each combination of s and q, we 
allow the angle of the source trajectory, a, to vary, seeding each run with values of a from to 360 



degrees in steps of 5 degrees. For our models, we use the map-making method of lDong et al.l (120061 ) 
when the source is within two source radii of the position of the center of magnification. Outside 
this time range, we use the hexadeca pole or quadrupole approximations for the magnification 



(jPeicha fc Hevrovskvll2009l : lGouldll2008h . We used a Markov Chain Monte Carlo to find the best-fit 



parameters and uncertainties for each s, q combination. 

The grid search reveals an overall improvement in x 2 relative to the point lens model. We find 
three % 2 minima for different angles for the source trajectory. For central caustics with planetary 
mass ratios, the caustic is roughly triangular in shape with a fourth cusp where the short side of the 
triangle intersects the binary axis; the three trajectories roughly correspond to the three major cusps 
of the caustic. An example caustic is shown in Figure H] along with the trajectories corresponding 
to the three minima. The angles of the three trajectories are approximately a = 0, 115, and 235 
degrees, and we will refer to them as trajectories "A", "B", and "C", respectively. 

We then refine our grid of s and q around each of these three minima. We repeat these 
fits accounting for various microlensing degeneracies. First, we fit without parallax and assuming 
s < 1. Then, we add parallax and fit both with uq > and uq < to see if this degeneracy is 
broken. Finally, we fit 2-body lens models with s > 1 and no parallax, since there is a well known 
microlensing degeneracy that takes s — > s _1 . 

The best-fit solution has % 2 = 6637.96 and a = 236.4. This reflects an improvement in % 2 of 
A X 2 ~ 140 over the point lens solution. There is no preference for s < 1 over s > 1, but trajectory 
C is preferred by Ax 2 > 35 over trajectories A and B. The parameters and their uncertainties for 
the planet fits are given in Table The mass ratio between the lens star and its companion is firmly 
in the planetary regime: q = 10~ 3 01 . Furthermore, planetary mass ratios are clearly preferred 
over "stellar" mass ratios (q ~ 0.1), which are disfavored by more than A% 2 = 60. Parallax further 
improves the fit by only Ax 2 ~ 10 and has little effect on the other parameters. 
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To compare the point lens and planetary models, in Figure [5j we plot the "x 2 residuals" : the 
difference between the cumulative x 2 distribution and the expected value Xj = i.e. each 
point is expected to contribute x] — 1- For both the point lens and the planet fits, the distribution 
rises gradually over the peak of the event. This is expected since 1-parameter errors do not account 
for correlated noise. However, in the x 2 residuals for the point lens, there is a jump seen at the 
time of the first limb-crossing. This jump is even more pronounced when looking at the difference 
between the planet and point lens models. The jump is caused by MOA data at the time of the 
limb-crossing that do not fit the point lens well, thereby causing an excess increase in x 2 ■ This is 
exactly the time when we expect to see planetary signals. 

4.4. Reliability of the Planetary Signal 

Although Ax 2 ~ 140 appears to be significant, we are hesitant to claim a detection of a planet. 
The planetary signal is more or less equally divided between the jump at the first limb-crossing 
and a more gradual rise after the second limb-crossing (see the third panel of Fig. [5] showing the 
difference between the point lens and planet models). One could argue that the gradual rise, due 
to a slight asymmetry in the planet light curve, could be influenced by large-scale correlations in 
the data. Comparing Figures [3] and shows that most of the signal at the first limb-crossing comes 
from only a few points. A careful examination of the residuals in Figure [3] shows that while the 
residuals to the planet model are smaller than for the point lens model, they are not zero, and 
the didactic residuals do not go neatly th rough the difference between the models as they do for 



MOA-2008-BLG-310 (jJanczak et al.ll2010i ). Hence, the evidence for the planet is not compelling. 



If we repeat the analysis using 2-parameter errors, we find a similar planetary solution, although 
the exact values of the parameters are slightly different The total signal from the planet is 
significantly degraded for 2-parameter errors, with only A% 2 ~ 80 between the best-fit planet and 
point lens models. Table [3] gives parameters for the complete set of point lens and planet fits for 
2-parameter errors. The residuals and error bars over peak may be compared to 1-parameter errors 
in Figure O 

We also show the x 2 residuals for 2-parameter errors in the bottom set of panels in Figure 
[5j They are more or less flat over the peak, showing that they track the data well in this region. 
The offset from zero is caused by systematics elsewhere in the light curve. The difference plot 
(bottom-most panel) shows that the planet fit is still an improvement over the point lens fit, but 
the signal from the planet at the first limb-crossing is much weaker. This is a natural consequence 
of 2-parameter errors, since the data at the peak of the event, where the planetary signal is seen, 



For 2-parameter errors, model "A, wide" appears to be competitive with model "C". However, this solution 
requires that the source trajectory pass over the planetary caustic at the exact time to compensate for a night for 
which the MOA baseline data are high by slightly more than la compared to other nights at baseline. If the data 
from this night are removed, the remaining data predict a different solution with the planetary caustic crossing 18 
days earlier. Because this solution is pathological, we do not consider it further. 
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have much larger renormalized error bars than for 1-parameter error 



Regardless of the error renormalization, this planetary signal is smaller than the Ax 2 of any 
securely detected high- magnification microlensing planet. Previously, the smallest Ay 2 ever re- 
porte d for a high-magni fication event was for MOA-2008-BLG-310 with A X 2 = 880 d.Tanczak et al 



20ld ). lYee et all (|2012l ) discuss MOA-2011-BLG-293, an event for which the authors argue the 
planet could have been detected from survey data alone with Ax 2 = 500. However, although the 
planet is clearly detectable at this level, it is unclear with what confidence the authors would have 
claimed the detection of the planet in the absence of the additional followu p data, whi c h incr eases 
the significance of the detection to Ax 2 = 5400. At an even lower level, iRhie et al.l (2000J) find 
that a planetary companion to the lens improves the fit to MACHO- 9 8-BL G-35 at Ax 2 = 20, but 
they do not claim a detection. As previously mentioned, iGould et al.1 (120101 ) suggest the minimum 
"detectable" planet will have 350 < Ax 2 < 700. However, this threshold has not been rigorously 
investigated; the minimum Ax 2 could be smaller. 

Because of the tenuous nature of the planetary signal, we do not claim to detect a planet in 
this event, but since including a planet in the fits significantly improves the X 2 , we will refer to this 
as a "candidate" planet. 

Finally, it is interesting to note that even though the Ax 2 for the planetary model is too small 
to be considered detectable, for both 1-parameter and 2-parameter errors the parameters of the 
planet (s and q) are well defined (see Tables [2] and [3]). Central caustics can be degenerate, especially 
when they are much smaller than the source size, so we might ex pect a wide range of possible mass 
ratios in this case since the limb-crossings are not well-resolved (|Hanll2009i ). However, perhaps we 
should not be surprised that the planet parameters are well-constrained: A^ 2 = 140 is formally a 
highly significant signal, which would plausibly lead to reasonable constraints on the parameters. 
In this case, because we believe that the signal could be caused by systematics, by the same token, 
the constraints on the parameters may be over-strong. We conjecture that the limb-crossing signal 
does not constrain q and that this constraint actually comes from the asymmetry of the light curve, 
since small, central caustics due to planets are asymmetric whereas those due to binaries are not. 

5. # E and /x rel 

Because the source size, p, is well measured, we can determine the size of the Einstein ring, #e 5 
and the relative proper motion between the source and the lens, p re \ from the following relations: 



P 



and 



0E 

Mrel — — ■ 
EE 



(5) 



Keeping the limb-darkening parameters fixed, we find the best fit for the normalized source size 
to be p = (2.70 ± 0.06) x 10~ 3 for the planetary fits; the value is comparable for the point lens 



4 Note that while the outliers are slightly different for 1-parameter and 2-parameter errors, no points were rejected 
in either case during the first limb- crossing, 5365.13 < t(HJD') < 5365.18, when the main planetary signal is observed. 
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fits. We convert the (V — I) color to (V — K) using iBessell &: Brettl ()1988l ) and obtain the surface 
brightness by adopting the relation derived by iKervella et al.l (|2004l ) . Combining the dereddened I 
magnitude with this surface brightness yields the angular source size 6* = 1.03±0.07 //as. The error 
on 6* combines the uncertainties from three sources: the uncertainty in flux (f s ,i), the uncertainty 
from c onverting (V — I) color to the surface brightness, and the uncertainty from the iNataf et al. 
(|2012l ) estimate of lo,ci- The uncertainty of f s j is 3%, which is obtained directly from the modeling 



output. We estimate the uncertainty from the other factors (Z) to be 7%. The fractional error in 
(/* is given by [(l/4)(a /s/ ) 2 + (az/Z) 2 } 1 / 2 = 7%, which is also the fractional error of the proper 
motion fi and the Einstein ring radius #e (|Yee et al.l 12009 ). Thus, we find #e = 0.38 ± 0.03 mas 
and /i re i = 7.1 ± 0.6 mas yr _ . 

6. The Lens as a Possible Member of NGC 6553 



This microlensing system is close in projection to the globular cluster NGC 6553. The cluster 
is at (R.A., decl.) = (18 h 09 m 17?60, -25°54 / 31. // 3) (J2000 .0), with a distance of 6.0 kpc from the 
Sun and 2.2 kpc from the Galactic Center (lHarrislll996l ). The half-light radius of NGC 6553 
is 1.03' (lHarrislll996l ). which puts the microlensing event 6.5 half-light radii (6.7') away from the 
cluster center. By plotting the density of excess stars over the background, we find that about 6% 
of the stars at this distance are in the cluster (Figure [6]) . 

Whether or not the lens star is a member o f the c luster can be constrained by calculating 
the proper motion of the lens star. IZoccali et al.1 (|200ll ) found a relative proper motion of NGC 
6553 with respect to the Bulge of /x = (/i,;, /Xf,) = (5.89,0.42) mas yr" 1 . The typical motion of 
Bulge stars is about 100 km s" 1 corresponding to about Smasy r -1 given an estim ated distance of 



7.7 kpc toward the source along this line of sight (I = 5.1°; INataf et al.ll2012l ). Therefore the 
expected amplitude of the lens-source relative proper motion if the lens is a cluster member is 
/^rei = 7 ± 3masyr _1 , which is consistent with the measured value in Section ([5]). 

Combining the measurement of the stellar density with the proper motion information, we 
find the probability that the lens is a cluster member is considerably higher than the nominal 
value based only on stellar density. First, of order half the stars in the field are behind the source, 
whereas the lens must be in front of the source. Second, the lens-source relative proper motion is 
consistent with what would be expected for a cluster lens at much better than 1 a, which is true 
for only about 2/3 of events seen toward the Bulge. Combining these two effects, we estimate a 
roughly 18% probability that the lens is a cluster member. 

One way to resolve this membership issue is by measuring the true proper motion of the lens 
as was done for a microlensing event in M2 2 for which the lens was confirmed to be a member of the 
globular cluster (|Pietrukowicz et al.ll2012l ). For this event, we have calculated the relative proper 
motion of the lens and the source to be 7.1 ± 0.6 mas yr -1 . This is consistent with the expected 
value if the lens were a member of the cluster. About 10 years after the event, the separation 
between the lens and the source star will be large enough to be measured with HST. Based on 
this followup observation, one will be able to clarify whether the lens is a member of the cluster by 
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measuring the vector proper motion. If it is a member of the cluster, its mass may be estimated 
from 

q2 

Mi ens = (6) 
K7T re l 

where k = 4G/c 2 AU = 8.14 mas/M© and 7r re i is the source-lens relative parallax. In addition, the 
metallicity of the lens could be inferred from the metallicity of the globular cluster. 

7. Discussion 

We have found a candidate planet signal in MOA-2010-BLG-311. The evidence in support of 
the planet is 

1. The planet substantially improves the fit to the data, 

2. In addition to a general improvement to the light curve, the planet produces a signal when 
we most expect it, i.e. the time of the first limb-crossing, 

3. The solution has a well-defined mass ratio and projected separation for the planet (excepting 
the well known s — > s _1 degeneracy). 

The magnitude of the signal depends on whether the error bars are renormalized using 1-parameter 
or 2-parameter error factors. We argue that for limb-crossing signals, which are very short duration, 
1-parameter errors are more appropriate since the signal is likely to be on a scale that is smaller than 
the correlation length of the correlated noise. With 1-parameter errors, we get a larger planetary 
signal, Ax 2 = 140, but one that is still too small to claim as a secure detection. 

Examining the residuals to the light curve and the x 2 residuals shows that the planetary signal 
is dispersed over many points at the peak of the light curve. It comes from an overall asymmetry 
near the peak plus a few points at the limb-crossing. Because the signal is the sum of multiple cases 
of low-amplitude deviations, it is plausible that the microlensing model could be fitting systematics 
in the data, which is why the planet signal is not reliable. 

This event suggests that central-caustic (high-magnification) events and plan etary-caustic 



events require different detection thresholds. The detection threshold suggested in I Gould et al 
(201CJ) of 350 < Ax 2 < 700 was made in the context of high-magnification events, and our expe- 



rience so far is consistent with th is. A plane t was clear l y det ectable in MOA-2008-BLG-310 with 
Ax 2 = 880 ( Janczak et al.ll2010l ). However, lYee et al.l ( 20121 ) are uncertain if a planet would be 



claimed with Ax 2 = 500 for MOA-2011-BLG-293. Here, Ax 2 = 140 is clearly insufficient to detect 
a planet. In contrast, the planeta ry caustic crossing in OGLE-2005-BLG-390 produced a signal of 



Ax 2 = 532 (jBeaulieu et al.1 120061 ) . and the planet would probably still be detectable if the error 
bars were 50% larger (Ax 2 ~ 200) and might even be considered reliable if the error bars were 
twice as large. 



We suspect the reason for the different detection thresholds is that the information about the 
microlens parameters and the planetary parameters comes from different parts of the light curve. 
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For planetary-caustic events, the planet signal is a perturbation separated from the main peak. 
Thus, the microlens parameters can be determined from the peak data independently from the 
planetary parameters, which are measured from the separate, planetary perturbation. In contrast, 
for high-magnification events, the planetary perturbation occurs at the peak of the event, so the 
microlens and planetary parameters must be determined from the same data. 

The detection threshold for planetary-caustic events will have to be investigated in more detail. 
If it is truly lower than for high- magnification events, this is good news for second generation 
microlensing surveys since that is how most planets will be found in such surveys. At the same 
time, it points to the continued need for followup data of high-magnification events since these 
seem to have a higher threshold for detection, requiring more data to confidently claim a planet. 
This is an important consideration because high-magnification events can yield much more detailed 
information about the planets. 



Additionally, Irian Kiml (|2009l ) show that the magnitude of the planetary signal should de- 
crease as the ratio between the caustic size and the source size (w/p) decreases for the same 
photometric precisioiJE There are several cases of events for which w/p < 1: in this case, we 
have w/p = 0.12 and Ay 2 = 140; th e brown dwarf in MOA-2009-BLG-411L has w/p = 0.3 

= 1070 
= 880 



Janczak et al 



and Ay 2 = 580 (jBachelet et alJ 12012): MOA-2007-BLG-400 has w/p = 0.4 and A X 2 



Dong et al.ll2q0flh: a nd in the case of MOA-2008-BLG-310, the value is w/p = 1.1 with Ax 2 



2010) . This sequence is imperfect, but the ph otometry in t he fou r cases is far from 



uniform, and it seems that in general the trend suggested by lHan &; Kiml (|2009l ) holds in practice. 



Finally, we note that a large fraction of the planet signal comes from the MOA data, but 
in the original, real-time MOA data, this signal would not have been detectable since the data 
were corrupted over the peak of the event. It is only after the data quality was improved by 
rereductions, which in turn were undertaken only because the event became the subject of a paper, 
that we recovered the planet candidate. This points to the importance of rereductions of the data 
when searching for small signals. In the current system of analyzing only the events with the most 
obvious signals, this is not much of a concern. However, if current or future microlensing surveys 
are systematically analyzed to find signals of all sizes, this will become important. 
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Fig. 1. — Color-magnitude diagram (top) and color-color diagram (bottom) of the field of view of 
MOA-2010-BLG-311 constructed from the CTIO observations. The small black dots represent the 
field stars; stars used in constructing the color-color diagram are shown as the green squares in 
both panels. The black line in the bottom plot shows the measured instrumental (I — H) color of 
the source. The (V — I) color of the source is then derived from the intersection of this line with 
the (V — I) — (I — H) relation established by the field stars (green squares), producing the black 
dot in the bottom panel. The position of the source in the CMD (large black dot) is given by this 
(V — I) measurement and the source flux from the fits to the light curve. The large red dot shows 
the centroid of the red clump. 
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Fig. 2. — Light curve of MOA-2010-BLG-311. Data from different observatories are plotted in 
different colors. The data from Bronberg {medium pink) have been binned for clarity in the figures; 
only unbinned data were used in the fitting. The black line shows the best-fit point lens model; 
on this scale, the best-fit planetary model appears very similar. The error bars reflect 1-parameter 
errors (see Sec. 12. 3|) . 
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Residuals: 1— parameter errors 
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Fig. 3. — Residuals over the peak to the best-fit point lens and planetary microlensing models. The 
top set of panels shows 1-parameter errors, and the bottom set shows 2-parameter errors. Note that 
for the time range shown, which is the same as for the bottom panel for Fig. [21 the error bars tend 
to be larger for 2-parameter errors than for 1-parameter errors. Careful inspection of the residuals 
to the point lens and planet models shows improvement around the times of the limb-crossings of 
the source star and also on the falling side of the light curve. In the middle panel of each set, the 
red lines show the difference between the best planetary and best point lens models. The didactic 
residuals in the middle panel are the sum of the red line and the residuals to the best planet model 
(bottom panel). The dashed lines indicate the limb-crossing times. The dotted lines in the top 
panels of each set indicate the time ranges shown in the bottom panels. Note that the scales for 
the top and bottom residuals panels are not the same. The colors of the points are the same as in 
Fig. [2j The Bronberg data have been binned for clarity (as in Fig. [2]) 
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Fig. 4. — The caustic (red) and three source trajectories corresponding to the three minima dis- 
cussed in Section 14.31 Trajectory C (black) is preferred by Ax 2 > 35 over trajectories A (green) 
and B (blue) for models with 1-parameter errors and without parallax. The source is shown to 
scale as the large circle; the lines with arrows indicate the trajectories of the center of the source. 
The abscissa in these plots is parallel to the binary axis with the lens star close to the origin and 
planet to the right. The scale is such that 1.0 equals the Einstein radius. 
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Fig. 5. — The difference between the cumulative x 2 distribution and the expected distribution 
(52f Xj = -^0) i- e - the x 2 residuals, for both 1-parameter (top) and 2-parameter (bottom) errors. 
The top two panels of each set show the distributions for the point lens model (a) and the planetary 
models (b). The bottom panels (c) of each set show the difference between the top and middle 
panels. The distributions for each data set are plotted separately and are shown over the same 
time range as the bottom panel of Fig. [2j data sets without points in this time range are not 
shown. The colors are as in Fig. [2l The thick black line shows the total distribution for all data. 
The limb-crossing times are indicated by the dashed lines. Note the jump in the MOA data (light 
green) at the time of the first limb-crossing (HJD' ~ 5361.5). The signal is much more pronounced 
for 1-parameter errors than for 2-parameter errors. Note that the vertical scales in the two sets of 
panels are different. 



Fig. 6. — DSS image of NGC 6553 (left). The position of the microlensing event is indicated by the 
circle. Excess star density over the background around the globular cluster NGC 6553 (right). The 
center of the cluster is placed at r = 0. The microlensing event is at r = 6.7' (dashed line). The 
dip at ~ 6' is caused by a quasi-circular dust lane, which may be seen in the image. Even though 
the density of stars drops quickly as a function of radius, cluster members still comprise about 6% 
of the stars at r = 6.7'. 
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Table 1. Data Properties for Two Methods of Error Renormalization. The observatories are 
listed in order of longitude starting with the most Eastward. If data were taken in more than one 
filter at a given site, different filters are given on successive lines. The error renormalization 
coefficients and method for removing outliers are described in Section! 



1-parameter Errors 2-parameter Errors 
Observatory Filter k N dSLtSL k e min iV data 



Mt. Bigelow 


I 


1.63 


44 


1.63 


0.0 


44 


Molehill 


Unfiltered 


0.72 


69 


0.72 


0.0 


69 


Kumeu 


Wratten #12 


1.19 


188 


1.18 


0.0 


188 


Farm Cove 


Unfiltered 


1.27 


52 


1.26 


0.0 


52 


Auckland 


Wratten #12 


0.98 


84 


1.00 


0.0 


84 


Vintage Lane 


Unfiltered 


4.87 


112 


3.05 


0.004 


112 


B&C 


I 


4.07 


132 


1.01 


0.025 


136 




V 


1.15 


53 


0.66 


0.03 


55 


MOA 


MOA-Red 


1.68 


4452 


1.55 


0.003 


4434 


Canopus 


I 


3.02 


28 


2.87 


0.0 


29 


Wise 


Unfiltered 


0.52 


70 


0.53 


0.0 


70 


Bronberg 


Unfiltered 


1.26 


727 


1.27 


0.0 


727 


SAAO 


I 


2.60 


128 


2.21 


0.0015 


127 


Liverpool 


SDSS-i 


1.85 


120 


1.04 


0.007 


119 


La Silla 


I 


10.02 


169 


3.79 


0.004 


174 


CTIO 


I 


1.33 


22 


1.34 


0.0 


22 




V 


0.50 


3 


0.50 


0.0 


3 








74 






74 


OGLE 


I 


1.36 


429 


1.32 


0.008 


429 



a These data were not used in the modeling. They were only used to determine the source color 
(see Sec. [3]). 



Table 2. Fits with 1-Parameter Errors. The A% 2 is given relative to the x 2 °f the best-fit 
planetary model with s < 1 and without parallax (x 2 = 6637.96), i.e. model "C"; positive 
numbers indicate a worse fit and negative numbers indicate an improvement relative to this 
model. The point lens models are given first, followed by the planetary models; "A", "B", and 
"C" denote the three planetary models with distinct values of a corresponding to the three x 2 

minima. 



Model 




Ax 2 


t - 5365. 


Uq 




P 


a 


log s 


log q 


1"B,JV 










(HJD') 




(days) 




(°) 










Point Lens 




136.44 


0.19615(4) 


0.00152(3) 


20.34(42) 


0.00260(5) 












PL, parallax 




69.61 


0.1978(2) 


0.00167(4) 


19.37(41) 


0.00275(6) 








3.16(41) 


-1.34(20) 


PL, parallax, 


-n 


123.89 


0.19613(9) 


-0.00158(4) 


20.23(42) 


0.00262(5) 








-1.07(36) 


-0.98(27) 


A 




35.26 


0.19615(4) 


0.00167(3) 


19.04(34) 


0.00279(5) 


347.7(6) 


-0.12(1) 


-4.46(8) 






A, parallax 




-7.02 


0.1976(2) 


0.00172(4) 


18.92(40) 


0.00283(6) 


347.4(3) 


-0.08(1) 


-4.9(1) 


2.82(44) 


-1.29(22) 


A, parallax, - 


-n 


-3.98 


0.19626(9) 


-0.00184(5) 


18.51(38) 


0.00289(6) 


-346.(1) 


-0.17(2) 


-4.17(8) 


-2.11(46) 


-2.43(40) 


A, wide 




33.07 


0.19615(4) 


0.00167(4) 


19.06(39) 


0.00279(6) 


348.1(6) 


0.12(1) 


-4.44(8) 






B 




49.06 


0.19614(4) 


0.00159(4) 


19.71(43) 


0.00269(6) 


118(1) 


-0.43(4) 


-3.5(1) 






B, parallax 




21.45 


0.1973(2) 


0.00166(4) 


19.40(41) 


0.00275(6) 


119(2) 


-0.40(5) 


-3.6(2) 


2.21(45) 


-1.05(21) 


B, parallax, - 


-Uq 


31.33 


0.19617(9) 


-0.00169(4) 


19.46(41) 


0.00274(6) 


-115(1) 


-0.40(3) 


-3.5(1) 


-1.45(41) 


-1.50(35) 


B, wide 




49.11 


0.19615(4) 


0.00159(3) 


19.76(38) 


0.00268(5) 


118(1) 


0.43(4) 


-3.5(1) 






C 




0.00 


0.19613(4) 


0.00159(3) 


19.68(41) 


0.00270(6) 


236.4(7) 


-0.26(4) 


-3.7(1) 






C, parallax 




-10.65 


0.1968(3) 


0.00164(4) 


19.34(39) 


0.00275(6) 


235(1) 


-0.4(1) 


-3.4(3) 


0.89(51) 


-0.78(24) 


C, parallax, - 


-no 


-12.63 


0.19630(9) 


-0.00171(4) 


19.20(38) 


0.00278(6) 


-232(1) 


-0.51(8) 


-3.0(2) 


-1.11(46) 


-1.55(41) 


C, wide 




0.00 


0.19613(4) 


0.00159(3) 


19.73(39) 


0.00269(5) 


236.4(7) 


0.26(4) 


-3.7(1) 







Table 3. Fits with 2-Parameter Errors. The Ax 2 is given relative to the x 2 °f the best-fit 
planetary model with s < 1 and without parallax (x 2 = 6751.93), i.e. model "C"; positive 
numbers indicate a worse fit and negative numbers indicate an improvement relative to this 
model. The point lens models are given first, followed by the planetary models; "A", "B", and 
"C" denote the three planetary models with distinct values of a corresponding to the three \ 2 

minima. 



Model 


Ax 2 


to - 5365. 


"0 




P 


a 


log s 


log (J 


1"B,JV 








(HJD') 




(days) 




(°) 










Point Lens 


81.12 


0.19618(5) 


0.00152(3) 


20.51(44) 


0.00259(6) 












PL, parallax 


25.43 


0.1978(2) 


0.00164(4) 


19.80(43) 


0.00270(6) 








3.20(43) 


-1.14(21) 


PL, parallax, — uq 


72.23 


0.1961(1) 


-0.00160(4) 


20.23(44) 


0.00263(6) 








-1.26(42) 


-0.96(32) 


A 


6.37 


0.19611(5) 


0.00160(4) 


19.82(43) 


0.00269(6) 


354(2) 


-0.49(5) 


-3.1(1) 






A, parallax 


-12.98 


0.1972(3) 


0.00168(5) 


19.25(46) 


0.00278(7) 


349(3) 


-0.3(2) 


-3.7(6) 


1.89(70) 


-1.16(25) 


A, parallax, —no 


-8.07 


0.1963(1) 


-0.00172(5) 


19.27(42) 


0.00278(6) 


-348(2) 


-0.46(7) 


-3.1(2) 


-0.98(52) 


-1.59(48) 


A, wide 1 


-1.59 


0.19614(5) 


0.00157(3) 


20.05(37) 


0.00266(5) 


356(3) 


0.539(4) 


-3.01(3) 






B 


13.25 


0.19613(5) 


0.00160(3) 


19.79(41) 


0.00270(6) 


114(2) 


-0.51(6) 


-3.1(2) 






B, parallax 


-6.42 


0.1972(2) 


0.00165(4) 


19.54(41) 


0.00274(6) 


112(3) 


-0.5(1) 


-3.1(3) 


1.97(48) 


-1.02(23) 


B, parallax, —uq 


4.22 


0.1962(1) 


-0.00167(5) 


19.54(41) 


0.00274(6) 


-109(3) 


-0.50(6) 


-3.1(2) 


-0.99(47) 


-1.24(43) 


B, wide 


13.10 


0.19613(5) 


0.00159(4) 


19.82(42) 


0.00269(6) 


113(2) 


0.51(6) 


-3.1(2) 






C 


0.00 


0.19614(5) 


0.00158(3) 


19.94(40) 


0.00268(5) 


234(1) 


-0.40(7) 


-3.3(2) 






C, parallax 


-14.90 


0.1970(3) 


0.00164(4) 


19.51(41) 


0.00274(6) 


231(2) 


-0.56(8) 


-2.9(2) 


1.26(54) 


-0.96(24) 


C, parallax, — «o 


-10.15 


0.1963(1) 


-0.00166(5) 


19.49(42) 


0.00274(6) 


-230(2) 


-0.51(7) 


-2.9(2) 


-0.44(51) 


-1.04(46) 


C, wide 


-0.02 


0.19613(5) 


0.00158(3) 


19.97(42) 


0.00267(6) 


234(1) 


0.40(7) 


-3.3(2) 







1 This solution and its parameters should be treated with caution, since it corresponds to a pathological geometry. See 
footnote \3\ 



